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Abstract 

The experiments at LHC have shown that the SUSY (exotic) bound states are hkely to form 
bound states in an entirely similar fashion as ordinary quarks form bound states, i.e., quarkonium. 
Also, the interaction between two squarks is due to gluon exchange which is found to be very 
similar to that interaction between two ordinary quarks. This motivates us to solve the Schrodinger 
equation with a strictly phenomenological static quark-antiquark potential: V{r) = — Ar~^ -\- K^/r + 
Vq using the shifted large A^-expansion method to calculate the low-lying spectrum of a heavy 
quark with anti-sbottom (c6, bb) and sbottom with anti-sbottom (bb) bound states with is set 
free. To have a full knowledge on spectrum, we also give the result for a heavier as well as for 
lighter sbottom masses. As a test for the reliability of these calculations, we fix the parameters of 
this potential by fitting the spin-triplet (n'^Si) and center-of-gravity / 7^ experimental spectrum 
of the ordinary heavy quarkonia cc, cb and bb to few MeV. Our results are compared with other 
models to gauge the reliability of these predictions and point out differences. 
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I. INTRODUCTION 



Weak-scale Supersymmetric Standard Model (SSM) is the leading candidate for physics be- 
yond the Standard Model (SM) [1,2]. Supersymmetry (SUSY) is built on a solid theoretical 
and mathematical foundation. It is also well-motivated as an elegant solution to the gauge 
hierarchy problem and has merits of gauge coupling unification, dynanical electroweak sym- 
metry and providing a legitimate candidate for dark matter. SUSY predicts the existence 
of a super partner called SUSY particles (sparticles) corresponding to each ordinary par- 
ticle of SM. These sparticles should be accessible at the exist and constructing colliders 
such as Tevatron and LHC. Over the past years, great effort has been made to search for 
such sparticles. So far, no direct signal for SUSY has been observed and some lower mass 
bounds have been estabhshed for sparticles. The experimental results at LEP [3-8,9,10] 
and Tevatron [11-14], squarks must be heavier than about 100 GeV. However, most exper- 
imental searches for sparticles are performed with model-dependent assumptions and rely 
on a large missing energy cut. A long-hved hght SUSY bottom quark (sbottom), 6, and its 
anti-sbottom, 6, with a mass close to mt, (~ 4.9 GeV), roughly half the T(15') mass, has 
not been excluded by experiments [15,16]. Hence, a light sbottom and its anti-sbottom, are 
not excluded so far partly because of the ALEPH collaboration indication [3-8] and partly 
because of interesting scenario to explain the excess of bb pair production in hadron collisions 
than theoretical prediction by a factor two. Some analyses [15] showed that if the light b 
is an appropriate admixture of left-handed and right-handed sbottom quark, its coupling 
to Z boson can be small enough to avoid LEP-I Z decay bounds. In addition, a scenario 
with light gluino and long-lived light sbottom with mass close to the bottom quark was 
proposed in [16] with which the excess of measured bb pair production in hadron collision 
over QCD theoretical prediction by a factor two is explained successfully (cf. [17-19]). The 
data about bb pair production in hadron collision given by CDF and DO can be explained 
well by QCD theoretical production: e.g., we can learn the details from the web address in 
[17-19]. The CLEO exclusion of a 6 with mass 3.5 to 4.5 GeV [20] can also be loosed even 
avoided, since their analysis depends on the assumption for semi-leptonic decays of the light 
sbottom. Moreover, since sbottom is a scalar, based on the spin freedom counting only, its 
pair production rate at collision will be smaller than the bottom quark by a factor four, so 
the sbottom samples must be rarer than those of bottom quark in experiments. 
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In contrary, it is interesting to point out that some experiments seemingly favor such a 
hght sbottom. The ALEPH collaboration has reported experimental hints for a light sbottom 
with a mass around 4 GeV and lifetime of 1 ps [21]. A recent analysis of old anomaly in the 
MARK-I data for cross section of e"'"e~— s>hadrons shows that the existence of such a light 
sbottom can bring the measured cross section into agreement with the theoretical prediction 
[22]. As mentioned by Berger et al. [16], a light-gluino analysis was done by Baer et al. [23] 
in which the gluino is assumed LSP. Cheung and Keung [24] modified the analysis [23] by 
letting the light gluino decay into h and h and study the possible constraint and implication 
at LSP. Therefore, the light gluino and light sbottom scenario will certainly give rise to 
other interesting signals, e.g., decay of Xh into the light sbottom [25], enhancement of tthb 
production at hadron colliders [26], decay of T into a pair of light sbottoms [27] and flavor- 
changing effects in radiative decays of B mesons [28]. 

The phenomenology of a very light sbottom has been studied by many authors recently 
[24-26,28-32]. If such a hght sbottom indeed exist, new meson-hke bound states is formed 
by a pair of the sbottom and anti-sbottom [hh) and fermion-like ones by an ordinary quark 
with anti-sbottom [qb) [e.g., heavy quark q = c,b) may be also formed. 

Some quarkonium binding systems like cc, bb and cb have been studied with encouraging 
success [33-35], in the framework of the potential model using a strictly phenomenological 
static heavy quark- antiquark potential belonging to the generality V{r) — —Ar""" -\- kt^ -\- 
Vq, {a — 1,13 — 1/2). Hence, the parameters of this potential are fixed by fitting the 
experimentally measured triplet S-states and the center-of-gravity (c.o.g.) non S-states of 
cc and bb spectra to their calculated levels and taking into consideration their hyperfine 
splittings in the framework of non-relativistic quarkonium model. It has been found that 
the potential description is flavor- independent, i.e., the same potential describes equally well 
the cc and the bb systems. Therefore, if we take the potential parameters obtained from the 
fitting of the mass spectra cc and the bb systems, we may predict the Schrodinger bound- 
state masses of the exotic states. We take the very same values of potential parameters 
used for the observed bound states of ordinary quarks to predict the unknown exotic squark 
bound states. It is well-known that the interaction between two squarks is due to the gluon 
exchange [36] which is found to be very similar to that interaction between two ordinary 
quarks, and an interaction due to Higgs particle exchange [37]. Thus, the gluon exchange 
interaction in squarkonium motivates us to use very similar parameter set of present potential 



3 



model as in quarkonium [33-35]. 

Over the past years, the experiments at LHC have shown that the exotic bound states are 
hkely to form bound states in an entirely similar fashion as ordinary quarks form bound state, 
i.e., quarkonium. The long-lived sbottom is not excluded by conventional searches and an 
analysis should be done to verify that there are no additional constraints on the allowed range 
of sbottom masses and lifetimes. In addition, it is well-known that the interaction between 
two squarks is due to gluon exchange which is found to be very similar to that interaction 
between two ordinary quarks, and an interaction due to Higgs particle exchange. Thus, the 
gluon exchange interaction in squarkonium motivates us to use very similar parameter set of 
present potential model as in quarkonium. The purpose is to calculate the spectroscopy of 
(Qb) Q = c,b and {bb) in terms of potential model with Coulomb plus square-root potential 
[33-35] in which the parameters are fixed by heavy quarkonia (cc) and {bb) with sbottom 
mass, is set as a free parameter. Furthermore, in order to have a full knowledge on such 
a spectrum, we also give the result for a heavier sbottom masses from 3.0 GeV to 60.0 GeV. 

The outline of this paper is as following: In section 2, we first review the analytic solution 
of the Schrodinger equation for non-equal mass case for the spin-triplet S'-states and c.o.g 
(center of gravity) P- and D-states. In section 3, we briefly present the squarkonium produc- 
tion through the leptonic decay width and through the decay. We present the observed 
states and reproduce the calculated spectrum of the cc, cb and bb spectra in section 4. Also, 
the corresponding triplet -S'-states and the c.o.g. 1 — 1,2 for an ordinary SM heavy quark 
with anti-sbottom (cb, bb) and sbottom with anti-sbottom bb meson-like binding system with 
sbottom mass, is set as a free parameter. The conclusions are also given in section 5. 

II. SPIN-AVERAGED BINDING MASS SPECTRUM 

We limit our discussion to the following generality of potentials [33,38-49]: 



where A and k are positive constants whereas Vq is taking any sign. These static quarkonium 
potentials are monotone nondecreasing, and concave functions satisfying the condition [41- 



V{r) 



Ar-'^ + KT^ + a, /3 > 



(1) 



49] 



V\r) > and V"{r) < 0. 



(2) 
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At least ten potentials of this generality, but with various values of the parameters, have 
been proposed in the hterature (see, for example, Ref . [33] and references therein) . Motyka 
and Zalewiski [34,35] have also explored the quality of fit in the region < a < 1.2, 

< /9 < 1.1 of the a, (3 plane reasonably well with coordinates a = 1, f3 = 0.5. Therefore, 
the five parameters {A, K,Vo,mc,mb) are fixed, in fitting the cc and the bb experimental 
triplet states in flavor-independent model, as the below values: 

y^^^ ^ _0^325250 ^ Q 7Qg3g^ _ 0.78891, (3) 
r 

with fltted quark masses: 

= 1.3959 GeV, = 4.8030 GeV, (4) 

where V{r), and are all in units of GeV. Notice that for the c quark, nic is roughly half 
the J/t/j mass and for b quark it is roughly half of the T(15') mass. Thus, the potential model 
(3) is convincing as it approaches to the perturbative QCD formula in the short-distance 
region and approached to the confining potential in the long-distance region. Consequently, 
in short-distance region, this potential involves the r^^ (Coulombic part) corresponding to 
one gluon exchange which is approaching to the perturbative QCD formula. The linear 
confinement part of the potential is ~ r, as in Cornell potential [50,51], is not seen. Such 
a linearly rising potential [50,51] is capable of confining quarks permanently and it can 
give rise to spectrum of particles containing light quarks in rough accord with experiment 
[52,53]. Probably the heavy quarkonia like bb, d> {B^ meson) and bb are too small to reach 
sufficiently far into the asymptotic region of linear confinement. On the other hand, the 
charmonium cc is too large to reach sufficiently far into the confining potential in the large 
distance particularly for excited states near and above the open flavor threshold. Perhaps a 
more flexible potential would exhibit the linear part, but one may be observing an effect of 
the expected screening of the interaction between the heavy quarks by the light sea quarks 
[54]. 

We choose the corresponding spin triplet c.o.g. states for the practical reasons that the 
masses of the spin singlet pseudoscalars for the bottomonium hb arc currently unknown or 
very poorly measured [55] and the two unknown charmonium cc states in the 3S and 4S 
multiplets are the 3^ So and 4^ So pseudoscalars. 

For a system of two composite particles, we shall consider the D-dimensional space 
Schrodinger equation for any spherically symmetric central potential in (3). Using ipiy) = 



Yi,m{^:4>)u{r)/r^^ it is straightforward to find the I ^ radial wave equation (in the 
usual units h — c— 1) : 



1 [k - {1 - a)][k - {3 - a) 



+ 



+ V{r)} u{r) = Er,,iu{r) , k^D + 2l-a, (5) 



where u{r) is the radial wave function, E^^i is the Schrodinger binding energy of meson and 
a is a proper shift. We follow the shifted 1/N or 1/k expansion method by defining 



drir- J m\Q 
m=0 V u / ^ 

and the binding energy expansion 



(6) 



CO -p(2-m) 

k 



n.l 



E 

m=0 



Q 



E„ 



(7) 



-1/2, 



where X = k (r/ro — 1), is an arbitrary point where the Taylor's expansions is being per- 
formed around. Following Refs. [33,38-49], we give the necessary expressions for calculating 
the binding energies to the third order: 

P 



Eo = V{ro) + 



E,^(5 

E2^P 



nr + -\uj 



_ (2-a) 
8/x 



;i -a)(3-a) ^ ^(1)' 



16H 



To 



(8) 

(9) 
(10) 

(11) 



where a^^^ and a^^-* are two useful expressions given by Imbo et al. [56-58] and also the 
parameter k is 

^=^8//rgy'(ro). (12) 

Hence, the total binding energy of the three-dimensional {D = 3) Schrodinger equation to 
the third order is 



En,i = Viro) + IroV'iro) + \ 



;i-a)(3-a) 

mix 



a 



(2) 



k 



+ 



k 



(13) 
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and the shift parameter is 



a = 2 - {2nr + 1) 



3 + 



roV"{ro) 
V'iro) 



1/2 



(14) 



The root, tq, in Eqs. (13)- (14) can be found through the relation: 



1 + 21+ {2nr + 1) 



3 + 



roV'iro) 
V'iro) 



1/2 



[8i^rlV'{ro)Y^\ n„/ = 0,l,2,3, 



(15) 



where the radial number rir — n — 1 with n = 1, 2, 3, • • • is the principal quantum number. 
Once To is determined through Eq. (15), hence finding the Schrodinger binding energy of any 
quarkonium system from Eq. (13) becomes relatively simple and straightforward. Finally, 
the corresponding ordinary or exotic bound states become 

M{qiqj)ni = mq, + nig. + 2En,i, 

M{qiqj)ni = mq. + rUg. + 2E„^i, (16) 

where m^^ and m^^, are the composite masses of the quark with antiquark and m^, and m~ 
squark with anti-squark Details of the model and the method of solution may be found in 
Refs. [33,38-48]. 

Now let us turn to the investigation of the spin-spin interaction. It is well-known that the 
system under study is a nonrelativistic, the treatment is based on the Schrodinger equation 
with a Hamiltonian [59,60] 

Ho = -\^ + V{t) + Vss, (17) 

where Vss is the spin-spin contact hyperfine interaction which is one of the spin-dependent 
terms predicted by one-gluon exchange (OGE) forces [59,60]. Recently, the spin-spin part, 
in momentum space {q — /j), was found to be [61-63] 



Vss{mi,m2,q) = 
with Wilson coefficient 



-9siQ) 



C3{q,m) 



3mim2' 



and d{q) — 



—j^ — C3{q,mi)c3{q,m2) - GN^diq) 



(18) 



1 / as{m2] 



8N^ 



-9/25 



Qs(g) 



-18/25' 



(19) 

where Nc is the number of colors, gs{q) is the running coupling constant [61-65]. The formula 
(19) improves upon the one-loop perturbative calculation in two important respects: (i) it 
is independent of /i and (ii) also includes the higher order logarithmic terms. 



If the coefficients are calculated at tree level; i.e., C3(//, m) = 1, = 0, the potential 
reduces to the Eichten-Feinberg result [66,67]. And if these coefficients are expanded to order 
as{fi) then reduced to a one-loop quarkonium spin-spin interaction in the nonrelativistic case 
[68-70] which is responsible for the hypcrfine splitting of the mass levels [71-84] 

^ = sS^'W - i) • f^"' 

adapted from the Breit-Fermi Hamiltonian. The number | substituted from the prod- 
uct of the spins corresponds to the recent assumption that the unperturbed nonrelativistic 
Hamiltonian gives the energy of the triplet states. Since for the states with orbital angular 
momentum L > the wave function vanishes at the origin, the shift affects only the S states. 
Thus, the only first order effect of this perturbation is to shift to the pseudoscalar ^5*0 states 
down in energy: 

with the wave function at the origin is calculated by using the expectation value of the 
potential derivative via [33,41-46,49,82-84] 

An application of the last formula needs the value of the wavefunction at the origin. This 
can be achieved by solving the Schrodinger equation with the nonrelativistic Hamiltonian 
and the coupling constant. In such an approach, the QCD strong coupling constant as{4:ii'^), 
on the renormalization point /i^ is not an independent parameter. It can be connected (in 



the MS renormalization scheme) through the two-loop relation [34,35,85-87] 

, o, 27rr7^ - 1 2a 

MS 

where Po — 11 — Like most other authors (see, for example, Refs. [33-35,72-75]), 
the strong coupling constant ^^(m^), is fitted to the experimental charmonium hyperfine 
splitting numbers AHF(lS,exp) w 116.5 ± 1.2 MeV and AHF(2S,exp) w 48.1 ± 4 MeV [41- 
46,55], which yields 

as{ml) = 0.254. (24) 

Knowing the coupling at the scale m^, we obtain the couplings at other scales as follows. 
The number of fiavours (n/) is put equal to three for A/i^ < nil interested in the 
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region Afj? < m^) , equal to four for ml > 4//^ > ml and equal to five for 4//^ > mf (we are 
not interested in the region 4//^ > mf). Then the value of q;s(4//^ = m^) from (23) is used 

to calculate A^~^^ and A^"'''*. Using the known value of A^~^' and Eq. (23) we find the 
value 

a^iml) = 0.200 and as(4/^c6) = 0-224, (25) 

and 

as{m^) = as{ml) and tt^(4/i^=) = as{4:f4b)- (26) 

" CO 

We follow these analysis in order to calculate the experimental binding masses of the heavy 
quarkonia cc, c6 and bb. At the end, it is worth to note that the sbottom 6 is a scalar, there 
is no spin-spin interaction (hyperfine splitting) for (Qb) and (bb), where Q — c,b. 



III. SQUARKONIUM PRODUCTION 

A. Production Through the Leptonic Decay 

The squarks might be discovered by detecting ordinary quark-antiquark bound states as 
resonances at LHC. (This is of course one of the main ways of studying charm and bottom 
quarks). The bound states of squarks are narrow resonances depend primarily on their 
leptonic widths (unless the squark itself has a very large width). Therefore, our aim will 
be to compute Fg. The leptonic decay widths of the heavy quarkonia and squarkonium 
are proportional to the squares of the wave functions at the origin. Therefore, they are 
significant only for the S states. To compute the decay rate for this process, we use non- 
relativistic bound state techniques since the squarks are expected to be very heavy objects 
(massive particles) and therefore their bound states are not relativistic systems. For the 
cc and bb quarkonium and bb sbottomonium systems, we shall consider the decays of the 
n^Si (vector, J^'^ — 1 ) states decay into charged lepton pairs, e.g. e+e~ pairs are usually 
calculated from the QCD corrected Van Royen-Weiskopf formula 

r.(.=5,^70 = We«(l-^). (27) 

where |'0(O)| is the bound state radial wavefunction at the origin. My the mass of the 
bound triplet (vector) state, a the fine-structure constant and Sg the charge of the quark 
in units of the electron charge. In the computation we have taken for ordinary quarks 
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ec = I and = — |, however, for squarks = |- For the cb quarkonium and cb,bb 
squarkonium, we consider the decays of the ti^Sq (pseudoscalar) states into ru^ pairs. Since 
the probabihty of such decays contains as a factor the square of the lepton mass, the decays 
into hghter leptons are much less probable [33-35] . For vector mesons containing light quarks 
(squarks) this formula leads to paradoxes (cf., e.g., Ref. [33] and references therein). For 
quarkonia, however, the main problem seems to be the QCD correction. Thus, in order to 
get quantitative predictions it is necessary to include higher order corrections which are not 
known. In order to estimate the missing terms we tried two simple forms. Exponentialization 
of the first correction 

C,(a,(mJ))=exp(-15^), (28) 

and Padeization 

C^M^D) = (29) 

37r 

We use the average of these two estimates as our estimate of the QCD correction factor 
extended to higher orders. The difference between Ci and C2 is our crude evaluation of the 
uncertainty of this estimate. Further, we have the relation 

where Cav is the averaged QCD correction factor. With our choice of parameters this formula 
reduces to 

_ 4771— 

r,(77^5i ^ II) = F(g)^A£;HF, (31) 

where niq (rriq) is the quark (squark) mass and F{c) — 7.07 x 10~^ and F{b) — F{b) — 
2 A3 X 10""^, see Eq. (34). The leptonic width is a small fraction of the total width. Experi- 
mentally, for narrow resonances (whose width is much smaller than the beam-energy spread) 
one measures the integrated area of the resonance cross-section. For a Breit-Wigner-type 

resonance this is connected to the leptonic width by the formula 



27r^(2J+l) rer, 
sdE- —, (32) 



and Ffe/r ~ 1 if the hadronic width predominates. The formula for the leptonic widths of 
the pseudoscalar cb quarkonium reads 

r... = ^fl \Vc,\' Ms^ml (1 - ^) (33) 
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where G is the Fermi constant, Vch ~ 0.04 is the element of the Cabibbo-Kobayashi- 
Masakawa matrix and the decay constant fs^ is given by the formula (cf., e.g., Ref. [33] and 
references therein) 



i2k'(o)r^2 



Ml 



C (aj, 



(34) 



where C{as) is QCD correction factor. Formally this decay constant is defined in terms of 
the element of the axial weak current 



A^(0) \B,{q)) = ifBjcbq^. 



(35) 



The QCD correction factor is 



C{as) = 1 



«s(/^cb) 



TT 



^ mb-rric, rub 

2 m — 

nib + fTic rric 



(36) 



With our parameters C{as) ~ 0.905 and since this is rather close to unity, we use it without 
trying to estimate the higher order terms. 
Let us note the convenient relation 



87ra,(4/i2J Mb, 
which for our values of the parameters yields 



(37) 



/s, = 65.2j^V^, 



(38) 



where all the parameters are in suitable powers of MeV. 



B. Production Through Decay 

We compute first the width of squarkonium J^'^ = 0++. The 0++ state decays almost 
entirely into two gluons and qq pairs. Explicit calculations give [89] 

r(o-^,9)=fi^^Vwo)r. (39) 



with m — ruq {rriq) is the constituent mass of quark (squark) of the bound state system and 



/5127ra,(m2)\ 47?^ ,,_^^^^2 
V J (1 + AR^y 



r(0- ^qq)^{ ^^^S^ 1 1^(0)1^ , (40) 
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with R — m/m, m — rUg {m-g) is the gluon (gluino) mass. In conventional quarkonium, 
one can extract the value of | '0(0)1^ from experimental data on the leptonic width of the 
quarkonium state (except, of course, for toponium). We choose to use a tentative input for 
|-(/'(0)|^ obtained by a Coulomb like potential (1/r gluonic behaviour), see Eq. (3). This is 
justified, since as the mass of the constituents of the bound state system goes higher (high 
mass regime), the short-range forces should approach the Coulomb like interaction whereas 
the confining hnear potential should be neghgible in the short distance hmit. We shall stick 
to a coulombic wavefunction, namely, [82-84] 



,2 1 / mas{rn?' 



\m\--[-''^^'q^] ■ (41) 



The relative branching ratio of Eqs. (39) and (40) is 



r(o++ 

Grifols and Mendez [89] computed r(0++ ^ gg) ^ 0.6 MeV and B > 3/8 for as(m|) 
0.1609, Trig 65 GeV and m = > 87 GeV or m < 12 GeV. 



IV. OBSERVED ORDINARY QUARK AND PREDICTED EXOTIC SQUARK 
SPECTROSCOPY 

The levels are labeled by S,P,D, corresponding to relative orbital angular momentum 
L = 0, 1, 2 between quark and antiquark. (No candidates for L > 3 states have been seen 
yet.) The spin of the quark and antiquark can couple to either S* = (spin singlet) or S* = 1 
(spin triplet) states. The parity of quark-antiquark or squark-anti-squark states with orbital 
angular momentum L is P = (—1)^+^, the charge-conjugation is C = (—1)^+'^. Thus, L — 
states can be ^Sq or ^5"!; L — 1 states can be ^Pi or ^-Po,i,2; L — 2 states can be or 
^-C^i,2,3- The radial quantum number is denoted by n [55]. 

The experimentally clear spectrum of relatively narrow states below the open-charm DD 
threshold of 3730 McV can be identified with the IS, IP and 2S cc levels predicted by 
potential models, which incorporates a color Coulomb term at short distances and a linear 
confining term at large distances [50,51]. 

A recent interest in charmonium spectroscopy [90] is revived because of the recent dis- 
covery of the long missing 77^(2^80) state of binding mass 3638 ± 4 MeV by the Belle Col- 
laboration [91,92], which has since then been confirmed by BABAR [93] and has also been 
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observed by CLEO in 77 collisions [94] . The observation of the 2^P2 state with binding mass 
3929 ± 5 MeV. The reported l^Pi he signal is in the decay chain t/j' — > 7r°/ic, he ^rjc- The 
masses found in two different inclusive analysis were 3524.8 ± 0.7 MeV and 3524.4 ± 0.9 
MeV (with an estimated systematic error of ~ 1 MeV) in exclusive decay (with six dif- 
ferent identified rjc final states.) Additional interest in cc spectroscopy has followed the 
discovery of the remarkable X(3872) by Belle [95] and CDF [96] in B decays to J/V^tt+tt" ; 
assuming that this is a real resonance rather than a threshold effect, the X(3872) is pre- 
sumably either a DD* charmed meson molecules [97-99] or a narrow J — 2 D-wme cc state 
[100,101]. Very recent observations of the X(3872) in 7J/V' and ajJ/t/j by Belle support 
l"*""*" DD* molecule assignment [102,103]. Experimental activity in the spin-singlet P-wave, 
with recent reports of the observation of exclusive l^Pj he state by CLEO measurement of 
mass 3524.65 ± 0.55 MeV [104,105]. The -0(3770) is generally assumed to be ^Di cc state, 
perhaps with a significant 2^Si component [106] . The four known cc states above the DD 
threshold, -0(3770), -0(4040), -0(4159) and -0(4415) are of special interest because they are 
easily produced at e+e~ machines. The mass -0(4040) which is a very interesting case for 
the study of strong decays. The 1 ?/'(4159) is 2^Di cc assignment. The final known above 
DD threshold is the 1 -0(4415) has the assignment 4^Si. 

On the other hand, the searched pseudoscalar bottomonium 7^(,(l^So) meson, the new 
observed meson T{4:^Si) = 10579.4 ± 1.2 MeV and the l^Dj (probably all or mostly J = 2) 
state with bound mass 10161.1 ± 1.7 MeV [55] and the pseudoscalar charmed bottom meson 
with an averaged experimental mass Bc{0'"^) — 6286 ± 5 MeV by Yao et al. have also 
revived this interest [55]. 

Wc first apply this model described above to the ordinary quarkonia known through the 
SM. Consequently, we calculate the cc, bb and cb{hc) quarkonium binding mass spectra in 
close agreement with above up-to-date experimental findings. The theoretically calculated 
quarkonium binding masses together with their S-states hyperfine splittings are listed in 
Tables 1-2. The hyperfine mass sphttings of the IS c6 state is predicted by the present 
potential model and other models are listed in Table 3. The calculated binding masses are 
found to be in close agreement to the recently observed ones. This allows us to extend this 
study to the unknown spectra of the squarks to predict their binding masses in a unified 
way. Further, in Tables 1-3, all hyperfine splitting calculations of the potential model try 
to reproduce the old experimental values, while the lattice calculations and perturbative 
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QCD favor the new values. No confirmed experimental data to check these predictions are 
available yet. The experimental finding Ahf(2S,cc) = ^'(S^Si) - 77^(2^80) w 48.093 ± 3.97 
MeV is found to be lower than the calculated values from the potential models [19,26]. In all 
cases, where comparison with the other models are significantly smaller than the splittings 
found by Eichten and Quigg [71] and Gupta and Johnson [107]. The QCD sum rules [108,109] 
finds the hyperfine splitting of the bottomonium AHF(lS,66,theory)= 6315^ MeV with the 
central value agrees well to several MeV with expectation. However, the uncertainty is 
too large to distinguish between the potential models. A lattice calculation [110] gives the 
hyperfine splitting AHF(lS,66,theory)= 60 MeV with a large uncertainty. The central value 
seems to be close to our model, but the uncertainty is big enough to be consistent with all 
the potential models quoted here. Furthermore, our model predicts nearly an approximate 
hyperfine splitting for the IS bottomonium and 2S charmonium as in the other potential 
models [34,35,49,71], lattice [111,112] and perturbation QCD [113,114]. 

The level fine and hyperfine splittings in charmonium and bottomonium together with the 
experimental and other models findings are listed in Table 4. The fine splitting in charmo- 
nium is found to be M^/(2S) — j^j/i/'(is)= 597 MeV. It is within 7.8 MeV from the experimental 
value. However, the splitting for bottomonium is found to be Mx'(2S) — Mr(is)= 571 MeV. 
It is within 8 MeV from the experimental value. 

Motivated by the great success of our earlier apphcations [33,38-49], we extend this study 
to produce the binding masses of the sbottom with anti-sbottom and heavy quark with 
anti-sbottom. In Table 5, we show the results about bb states. The numerical results about 
(qb) ,{q = b, c) states are shown in Table 6 and Table 7, respectively. Because sbottom is 
a spin zero particle, the spectrum of the corresponding bound states is simpler than that 
of {qq) states. We can see from these tables that when m^; » m^, the binding masses of 
the lowest states (bb) , (bb) and bb are very close to each other. This is reasonable because 
the strong interaction for these states is similar. Also, the leptonic decay widths of the 
sbottomonium system resonances for different values of the resonance sbottom mass and 
Cq = I are shown in Table 8. 

In this analysis, taking the unknown mass of the sbottom quark close to the ordinary 
bottom quark in the bb pair with = 4.5 MeV, wc find out the fine splittings for the 
S-states as follows: 569, 335, 252 MeV and for = 5.0 MeV as: 571, 332, 249 MeV. Thus, 
we remark that the fine splittings of the ordinary cc, bb and c6(&c) quarkonium systems and 
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the (qb), {q — b, c) and bb exotic quarkonium states are nearly same [130]. Finally, in general, 
the potential models seem to reproduce the experimental values much better. This feature 
would be understandable, since the potential models contain much more input parameters 
than the lattice or perturbative QCD models. 

We point out that Chang et al. [78] used a relativistic model to calculate part of these 
bound states for some J^^ — 0++, 1 ,2++, 3 corresponding to n^Po,n^Di,n^P2 and 
n^Dzi n — 0,1,2,3, respectively, which is entirely different than our c.o.g. calculations. 
Calculating states like n^Pi,n^P2 and n^D2 could help us to compare with other models 
[130]. 



V. CONCLUSIONS 



We have obtained the bound state masses and hyperfine energy splittings of a flavor- 
independent static quarkonium potential model for few ordinary quarkonium and scalar 
squarkonium mesons within the framework of the shifted A'"-expansion technique for L — 
0, 1, 2 states (see Tables 1-7). In Tables 1-4, we have shown the effectiveness of the employed 
static quarkonium potential model in producing the quarkonium bound-state masses to sev- 
eral MeV. Encouraged by this success of a flavor-independent potential model, we have also 
predicted the bound state masses of few unknown squarkonium systems for low- to- high 
squark masses (3.0 GeV-150.0 GeV) as shown in Tables 5-7. In finding the unknown squarko- 
nium energy splittings, we have used the quarkonium strong coupling constant as{m^) values 
to predict the squarkonium energy splittings. Essentially, this is because the type of interac- 
tion between two squarks is very similar to the interaction between two quarks (for reviews 
see, for example, [36] and references therein). Such an interaction is part of the one gluon ex- 
change interaction and is responsible for the mass differences. Apparently, the squarkonium 
fine and hyperfine splittings are found to be nearly same as their quarkonium counterparts 
if the squark mass is chosen near the quark mass (i.e., rriq ~ niq). Our conclusions are also 
consistent with the conclusions made by Ref. [130]. This is reasonable because the strong 
interaction for these states is similar (sec, for example, [130]). The exotic bound states are 
likely to form bound states in an entirely similar fashion as ordinary quarks form bound 
state, i.e., quarkonium [89]. Since the same potential model is used for conventional QCD 
bound states and sbottomonium, so it is obvious, from the present work and Ref. [130], 
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that the mechanism responsible for the binding of a bb couple is the same one responsible 
for the binding a scalar bb. In general, the type of interaction in squarkonium is very similar 
to that in quarkonium [36] . 

Furthermore, the calculation of the leptonic decay constant is important in predicting 
the cross-section. Therefore, we have calculated the leptonic decay constants in Table 8. We 
conclude that sbottomonium bound states can be detected as resonances at LHC. Above 
quark masses of about 2 GeV, Fg is much higher than the experimental upper limits even 
for squarks of charge |. Moreover, Fg increases rapidly with the squark mass. In addition, 
in the super symmetric front, squarks, sleptons and gauginos do also have a probability to 
be pairwise produced at LHC [89]. Squarks are likely to form bound states in an entirely 
similar fashion as conventional quarks from bound state, i.e., quarkonium [89]. The model- 
dependent squark leptonic decay widths are bigger than the ordinary quarks for < 
and > 10 GeV and give decay smaller than the energy splitting between bound states 
which is assumed to be a narrow resonance. As an illustrative example, with given squark 
mass rriq = 65 GeV and strong coupling constant as{m~) = 0.1609 in Ref. [89], we solve 
the Schrodinger equation for the strictly phenomcnological scalar potential (3) to obtain the 
squarkonium binding masses M{l^Si) = 127.745 GeV and M^I-^Sq) = 127.577 GeV for the 
singlet and triplet ground states, respectively. We, further, calculate the hyperfine splitting 
energy AE — 167.72 MeV and decay constant Fe = 0.556 KeV which is smaller than the 
energy-beam width. Hence, we conclude that this decay width is a wide resonance. On the 
other hand, the hidden supersymmetry showing up the squarkonium production through the 
Z° decay into a photon and the lowest lying state of the squarkonium system, the J^'-' = O^"*" 
state (^5*0 in spectroscopic notation). The integrated area of the resonance cross-section, 
Eq. (32), has a large (small) value for light (heavy) squark mass, respectively. For example, 
from Table 8, we find 0.008 x lO-^GeV'^ < f GresdE < 2.56 x 10"^ GeV-^ for the range 3.0 
GeV <m^< 150.0 GeV sbottom quark mass resonance. Hence, the detection of a low-mass 
squarks at LHC are much favored. 

Meanwhile, from the width of squarkonium 0++, Eqs. (39)-(42), we see that the relative 
branching ratio, B = 0.662 for = 60 GeV and m-g = 87 GeV by using the Coulomb 
potential result for |'?/^(0)|^ we get F(0+"'" — > 'gg) = 2.54 MeV with the strong coupling 
constant is taken relatively larger than Ref. [89]. For low- mass case — 10 GeV and 
rrig = 12 GeV, we obtain B = 0.498 4/8 which is above 3/8 and F(0++ ^ gg) ^ 0.424 
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MeV which is consistent with [89]. This result also favors low bottom squark mass [16]. 
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TABLE I: Observed and ealeulated binding mass spectrum of cc states (in MeV). AX denotes the 
mass shift of the spin-singlet state (n^So) from the spin-triplet state (n'^Si). 



State 


1 tPC\ 


Our work NR [59] GI [60] EQ [71] GJ [107] MZ [34,35] PDG [55] 




1(1--) 


3097 


3090 


3098 


3097 


3097 


3097 


3096 916 ± Oil 




1(0-+) 


—117 


— 108 


— 123 


—117 


— 117 


—117 


—116 516 ± 1 189 


IP (c.o.g) 


1(0,1,2)++ 


3521 


3594 3 


3525 


34Q9 


3526 


3521 


•3595 3-1-0 11 




2(1—) 


3694 


3672 


3676 


3686 


3685 


3690 


3686 093 ± 034 


A2^So 


2(0 +) 


-65.9 


-42 


-53 


-78 


-68 


-72 


-48 093 lb 3 966 


ID (c.o.g) " 


1(1,2,3)— 


3806 


3801 


3842 








3771.1 ± 2.4 


2P (c.o.g) " 


2(0,1,2)++ 


3944 


OA A O 

3943 


3963.3 








3929 ± 5 


V^(3-^5i) 


3(1-) 


4078 


4072 


4100 








4040±10 


2D (c.o.g) 


2(1,2,3)-- 


4150 


4161.2 


4211.4 








4159±20 


3P (c.o.g) 


3(0,1,2)++ 


4265 


4288.9 


4325.3 










V'(4'Si) 


4(1-) 


4377 


4406 


4450 








4415 ± 6 


SD (c.o.g) 


3(1,2,3)— 


4430 
















5(1-) 


4628 















n'^Di state. 
^2'^P2 state. 
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TABLE II: Observed and calculated binding mass spectrum of bb states (in MeV). AX denotes 
the mass shift of the spin-singlet state (n^So) from the spin-triplet state (n^Si). 



State 




Our work EQ [71] KR [115] MZ [34,35] PDG [55] 




1(1—) 


y4uu 


QA«A 




y4i:uu 


QAart on -I- n or 

y41:DU.oU ± U./D 


Al^So 


1(0-+) 




— o / 




f^R 7 


[LOU ) I 


IP(c.o.g) 


1(0,1,2)++ 


yyuu 


yo / O 


yyuo 


QQnn 
yyuu 


QOQQ 07 -1- n /I 1 

yoyy.o( ± u.4i 




2(1~) 


lUUol 






1 nn9'^ 


1 nn9'^ 9fi -1- n ^1 


A2^So 


2(0-+) 




AA 




98 




lD{c.o.g) {l^Dj 


state)** 1(1,2,3)— 


lUlOO 


1 ni 97 


lUlOD 


lUlOO 


1 m fil 1 -1-1 7 

iUiOi.i It 1. / 


2P(c.o.g) 


2(0,1,2)++ 


10261 


10231 


10259 


10260 


10260.237 ± 0.56 




3(1-) 


10364 


10339 




10355 


10355.2 ±0.5 


2L>(c.o.g) 


2(1,2,3)— 


10438 




10441 


10438 




3P(c.o.g) 


3(0,1,2)++ 


10527 




10520 


10525 






4(1-) 


10614 








10579.4 ± 1.2 


3L>(c.o.g) 


3(1,2,3)— 


10666 












5(1-) 


10820 











"Probably all or mostly J = 2. 
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TABLE III: The calculated be binding masses of the lowest S-states and its splitting compared 
with the other authors (in MeV). 

Work MB^il^SoY Mbj(135i) A15 

PDG (Yao et al.) [55] 6286 ± 5 - 

Our work 6290.8 6349 58.2 

Motyka and Zalewiski [34,35] 6291 6349 58 

Eichten and Quigg [71] 6264 6337 73 

Colangelo and Fazio [117] 6280 6350 70 

Chabab [109] 6250 ± 200 - 

Baker et al. [118] 6287 6372 85 

Roncagha et al. [119] 6255 6320 65 

Godfrey et al. [120] 6270 6340 70 

Bagan et al. [121] 6255 ± 20 6330 ± 20 75 

Brambilla et al. [122] - 632619^ 60 

Baldicchi and Prosperi [85-87] 6194 ~ 6292 6284 ~ 6357 65 < A15 < 90 

SLET [49]^* 62531^^ 63281^ 68 < A15 < 83 

SLET [49]= 6258i^i 63331^4 

Chen and Kuang [61,62] 6310 6355 45 

Gershtein et al. [123] 6253 6317 64 

Gupta and Johnson [107] 6267 6308 41 

"The averaged observed mass. 

^Averaging over the five values in Table 1 of [41]. 

'^We treat results of [71] in the same manner. 
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TABLE IV: Level hyperfine and fine splittings in charmonium and bottomonium (in MeV). 
Level splitting Our work [71]'^ [34,35]'^ [124-127]'^ [128]'^ [111]'' [129]^ [113,114]^= PDG [55] 





= M/;'(25) 




66 


78 


72 


98 


92 


43 




38 


48.093 lb 3.966 




= -^T(15) ■ 




58 


87 


57 


60 


45 




51 


44 


(160)? 


a2?(25) 


= Mx/(25) 




23 


44 


28 


30 


28 






21 







597 


589.177 ±0.023 


aP = M^,^:,s) - M^'(2S) 


384 


353.907 ± 9.966 




571 


562.96 ± 0.05 


A^p^^ = Mx'(35) - Mx'(25) 


333 


331.94 ±0.19 


A^**^ = -^T'(4S) - Afx'(3S) 


250 


224.2 ± 0.7 



"Potential model. 
''Lattice. 

^Perturbative QCD. 



TABLE V: The calculated binding energy masses of bb pair (in MeV). 



State mf' 


nS 


nP* 


nD" 




nS 


nP 


nD 


n 




1 3.0 


6033 


6455 


6718 


3.5 


6976 


7402 


7661 


n 




2 


6601 


6835 


7021 




7543 


7775 


7958 


n 




3 


6950 


7119 


7267 




7886 


8053 


8198 


n 




4 


7216 


7352 


7477 




8146 


8280 


8402 


n 




1 4.0 


7925 


8356 


8613 


4.5 


8880 


9316 


9572 


n 




2 


8493 


8724 


8904 




9449 


9679 


9858 


n 




3 


8831 


8997 


9139 




9784 


9948 


10088 


n 




4 


9087 


9219 


9339 




10036 


10167 


10284 


n 




1 5.0 


9839 


10280 


10535 5.5 


10800 


11248 


11501 


n 




2 


10410 


10640 


10817 




11375 


11604 


11780 


n 




3 


10742 


10905 


11043 




11705 


11866 


12003 


n 




4 


10991 


11120 


11236 




11951 


12079 


12193 


n 




1 6.0 


11764 


12219 


12471 8.0 


15639 


16121 


16373 


n 




2 


12343 


12572 


12747 




16236 


16468 


16639 


n 




3 


12671 


12831 


12967 




16560 


16719 


16851 


n 




4 


12915 


13041 


13155 




16797 


16921 


17031 


n 




1 10.0 19533 


20044 


20298 20.0 39120 


39793 


40066 


n 




2 


20152 


20388 


20558 




39879 


40142 


40314 


n 




3 


20476 


20634 


20764 




40215 


40379 


40506 


n 




4 


20709 


20831 


20938 




40442 


40563 


40666 


n 




1 40.0 78472 


79496 


79821 60.0 117887 119278 119663 


n 




2 


79559 


79884 


80072 




119329 119717 119925 


n 




3 


79942 


80126 


80259 




119766 119973 120114 


n 




4 


80179 


80307 


80410 




120019 120157 120264 



"The mass is in GeV. 

''The calculated c.o.g. binding energy mass for states n(0)++,n(l)++ and n(2) 
■^The calculated c.o.g. binding energy mass for states n(l) ,^^(2) and n(3) 
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TABLE VI: The calculated binding energy masses of cb pair (in MeV). 



State m^** 




nP 


nD 




nS 


nP 


nD 


n = 1 4.0 


5562 


5982 


6254 


4.5 


6052 


6471 


6743 


n = 2 


6140 


6381 


6576 




6629 


6869 


7063 


n = 3 


6504 


6682 


6837 




6992 


7168 


7323 


n = 4 


6785 


6929 


7061 




7271 


7415 


7546 


n = 1 4.8 


6346 


6766 


7037 


4.9 


6445 


6864 


7135 


n = 2 


6923 


7162 


7356 




7021 


7261 


7454 


n = 3 


7285 


7461 


7616 




7383 


7559 


7714 


n = 4 


7564 


7707 


7838 




7662 


7805 


7935 


n = 1 5.0 


6543 


6963 


7233 


5.5 


7036 


7456 


7726 


n = 2 


7119 


7359 


7552 




7611 


7851 


8044 


n = 3 


7481 


7657 


7812 




7972 


8148 


8302 


n = 4 


7760 


7903 


8033 




8250 


8393 


8523 


n = 1 6.0 


7530 


7950 


8219 


8.0 


9513 


9933 


10201 


n = 2 


8105 


8344 


8536 




10086 


10324 


10515 


n = 3 


8465 


8640 


8794 




10444 


10618 


10771 


n = 4 


8742 


8884 


9014 




10719 


10860 


10989 


n = 1 10.0 11502 11922 


12189 20.0 21479 


21900 


22166 


n = 2 


12074 12312 


12503 




22050 


22286 


22476 


n = 3 


12431 


12605 


12757 




22404 


22577 


22727 


n = 4 


12705 12846 


12974 




22676 


22815 


22942 


n = 1 40.0 41467 41888 


42153 60.0 61463 


61884 


62149 


n = 2 


42037 42273 


42462 




62033 


62269 


62457 


n = 3 


42390 42562 


42712 




62385 


62557 


62706 


n = 4 


42660 42799 


42925 




62655 


62794 


62920 



"The mass is in GeV. 
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TABLE VII: The calculated binding energy masses of bb pair (in MeV). 



State rnf' 


nS 


nP 


nD 


nS 


nP 


nD 


n = 1 4.0 


8695 


9130 


9385 4, 


.5 


9170 


9608 


9863 


n = 2 


9263 


9494 


9673 




9740 


9970 


10148 


n = 3 


9599 


9763 


9904 




10074 


10238 


10377 


n = 4 


9852 


9983 


10101 




10325 


10455 


10572 


n = 1 4.8 


9458 


9897 


10152 4, 


.9 


9554 


9994 


10248 


n = 2 


10028 


10258 


10435 




10124 


10354 


10531 


n = 3 


10361 


10524 


10663 




10457 


10620 


10759 


n = 4 


10611 


10741 


10857 




10707 


10836 


10953 


n = 1 5.0 


9650 


10090 


10345 5, 


.5 


10131 


10575 


10829 


n = 2 


10221 


10450 


10627 




10704 


10933 


11110 


n = 3 


10553 


10716 


10855 




11035 


11197 


11336 


n = 4 


10803 


10932 


11048 




11283 


11412 


11528 


n = 1 6.0 


10616 


11061 


11315 8, 


,0 


12567 


13022 


13274 


n = 2 


11189 


11419 


11594 




13145 


13375 


13550 


n = 3 


11520 


11681 


11819 




13474 


13634 


13770 


n = 4 


11767 


11895 


12010 




13717 


13844 


13957 


n = 1 10.0 14534 


14995 


15247 20.0 24457 


24935 


25187 


n = 2 


15117 


15347 


15520 




25051 


25283 


25454 


n = 3 


15444 


15603 


15738 




25376 


25534 


25667 


n = 4 


15685 


15811 


15923 




25613 


25737 


25847 


n = 1 40.0 44410 


44900 


45153 60.0 64393 


64888 


65140 


n = 2 


45013 


45246 


45417 




64999 


65233 


65404 


n = 3 


45337 


45495 


45627 




65323 


65481 


65613 


n = 4 


45572 


45696 


45805 




65558 


65681 


65789 



"The mass is in GeV . 
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TABLE VIII: Leptonic widths of the sbottomonium system resonances for different values of the 

resonance sbottom mass and e„ = i. 



State m{GeV) Te{KeV) m{GeV) Te{KeV) m{GeV) FeiKeV) m{GeV) FeiKeV) 



n = 


1 2.0 


1.715 


2.5 


1.635 


3.0 


1.571 


3.5 


1.523 


n = 


2 


0.694 




0.657 




0.620 




0.587 


n = 


3 


0.464 




0.445 




0.423 




0.401 


n = 


1 4.0 


1.487 


4.5 


1.461 


4.8 


1.450 


4.9 


1.447 


n = 


2 


0.557 




0.530 




0.516 




0.512 


n = 


3 


0.381 




0.362 




0.352 




0.349 


n = 


1 5.0 


1.444 


5.5 


1.433 


10.0 


1.504 


20.0 


2.032 


n = 


2 


0.507 




0.487 




0.381 




0.315 


n = 


3 


0.345 




0.330 




0.244 




0.174 


n = 


1 40.0 


3.405 


60.0 


4.869 


100.0 


7.863 


150.0 


11.641 


n = 


2 


0.324 




0.380 




0.533 




0.749 


n = 


3 


0.138 




0.131 




0.142 




0.172 
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